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Abstract 

In this paper we present further studies of recurrent configurations of Chip-firing games on Eulerian 
directed graphs (simple digraphs), a class on the way from undirected graphs to general directed graphs. 
A computational problem that arises naturally from this model is to find the minimum number of chips of 
a recurrent configuration, which we call the minimum recurrent configuration (MINREC) problem. We 
point out a close relationship between MINREC and the minimum feedback arc set (MINFAS) problem 
on Eulerian directed graphs, and prove that both problems are NP-hard. 

Keywords. Chip-firing game, critical configuration, recurrent configuration, Eulerian digraph, feed- 
back arc set, complexity, Sandpile model. 

1 Introduction 

A feedback arc set of a directed graph (digraph) G is a subset A of arcs of G such that removing A from 
G leaves an acyclic graph. The minimum feedback arc set (MINFAS) problem is a classical combinatorial 
optimization on graphs in which one tries to minimize \A\. This problem has a long history and its decision 
version was one of Richard M. Karp's 21 NP-complete problems (Kar72). The problem is known to be 
still NP-hard for many smaller classes of digraphs such as tournaments, bipartite tournaments, and Eulerian 
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multi-digraphs (CTY07; Hill; GHM07). We will prove that it is also NP-hard on Eulerian digraphs, a class 
in-between undirected and digraphs, in which the in-degree and the out-degree of each vertex are equal. 

Chip-firing game is a discrete dynamical system that has received a great attention in recent years, with 
many variants. The model is a kind of diffusion process on graphs that can be defined informally as fol- 
lows. Each vertex of a graph has a number of chips and it can give one chip to each of its out-neighbors 
if it has as many chips as its out-degree. A distribution of chips on the vertices of the graph is called a 
configuration. The model has been firstly defined on undirected graphs by A. Bjorner, L. Lovasz, and W. 
Shor (BLS91), and later on digraphs (BL92). The most important property of Chip-firing games is that if 
the game converges, it always converges to a unique stable configuration. This property leads to some re- 
search directions. A natural direction is the classification of all lattices generated by the converging games 
(LP01; Mag03). Most recently, the authors of (PP13) gave a criterion that provides an algorithm for deciding 
that class of lattices. In this paper we pay attention to another important direction initiated in a paper of N. 
Biggs. The author defined a variant of Chip-firing game on undirected graphs, the Dollar game (Big99), and 
studied some special configurations called critical configurations. Such configurations have many relations 
to algebraic and combinatorial objects such as spanning trees, acyclic orientations,Tutte polynomial, etc. A 
generalization to the case of digraphs was given by L. Levin et al. (HLMPPW08), where the authors defined 
recurrent configurations and presented many properties that are similar to those of critical configurations on 
undirected graphs. The authors also studied the Chip-firing game on Eulerian digraphs and presented several 
typical properties that can also be considered as natural generalizations of the undirected case. In this paper 
we continue this work and present generalizations of more surprising properties. 

A typical property of recurrent configurations is that any stable configuration being component-wise 
greater than a recurrent configuration is also a recurrent configuration. If the set of minimal recurrent con- 
figurations are known, one knows the set of all recurrent configurations. Hence it is worth studying properties 
of such recurrent configurations. It turns out from the study in (SchlO) that we can associate a minimal re- 
current configuration of an undirected graph G with an acyclic orientation of G. The acyclic orientations of 
G have the same number of arcs, namely \E(G)\, so do the total number of chips of minimal recurrent con- 
figurations. A direct consequence of this fact is that we can compute the minimum total number of chips of 
a recurrent configuration in polynomial time since we can compute easily a minimal recurrent configuration. 
It is natural to ask whether this problem can be solved in polynomial time for the case of digraphs. We will 
see that the problem becomes much harder than in the undirected case, even when the game is restricted to 
Eulerian digraphs with a sink. By giving the notion of maximal acyclic arc sets that can be regarded as a 
generalization of acyclic orientations of undirected graphs, we generalize the definitions and the results in 
(SchlO) to the class of Eulerian digraphs. Although natural, these generalizations are not easy to see from 
the studies on undirected graphs. They allow us to derive a number of interesting properties of feedback arc 
sets and recurrent configurations of the Chip-firing game on Eulerian digraphs, and provide a polynomial 
reduction from the MINREC problem to the MINFAS problem on Eulerian digraphs. We extend a result 
of (Flill) and show that the MINFAS problem on Eulerian digraphs is also NP-hard, which implies the 
NP-hardness of the MINREC problem on general digraphs. 

The paper is divided into two main sections. The first is devoted to the study of properties of the maximal 
acyclic arc sets that are complements of the feedback arc sets of an Eulerian digraph. The main result of 
this section is that finding an acyclic arc set of maximum size can be restricted to looking within particular 
subsets of acyclic arc sets. By using this result we prove that the MINFAS problem on Eulerian digraphs is 
NP-hard. It also gives a connection between the MINFAS problem and the MINREC problem on Eulerian 
digraphs that is presented in the second section. A direct consequence of this connection is the NP-hardness 
of the MINREC problem on general digraphs. 
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2 Acyclic arc sets on Eulerian digraphs 



Throughout this paper a graph always means a simple connected digraph. All results in this paper can be 
generalized easily to the case of multi-graphs. Traditionally, the vertex set and the edge set of a graph G are 
denoted by V(G) and E(G), respectively. An Eulerian digraph is a digraph in which the in-degree and the 
out-degree of each vertex are equal. An undirected graph is considered as a digraph in which for any edge 
linking u and v, we consider two arcs: one from u to v and another from v to u. With this convention an 
undirected graph is an Eulerian digraph. 

Let G = (V,E) be a digraph. For a subset A of £ let G[A] denote the graph (V',E') with V = V and 
E' — A. A feedback arc set F of G is a subset of E such that removing the arcs in F from G leaves an acyclic 
graph. An acyclic arc set A of G is a subset of E such that the graph G[A] is acyclic. Clearly, an acyclic 
arc set is the complement of a feedback arc set. A feedback arc set (resp. acyclic arc set) is minimum (resp. 
maximum) if it has minimum (resp. maximum) number of arcs over all feedback arc sets (resp. acyclic arc 
sets) of G. A feedback arc set A (resp. acyclic arc set A) is minimal (resp. maximal) if for any e e A (resp. 
e 6 E\A) we have A\{e] (resp. A U {e}) is not a feedback arc set (resp. acyclic arc set). 

From now until Proposition 1 we work with an Eulerian connected digraph G = ( V, E) (note that a 
connected Eulerian digraph is also strongly connected). A lot of properties of the acyclic arc sets of G are 
given in this section. The most important result is that finding a maximum acyclic arc set can be restricted to 
finding an acyclic arc set of the maximum size that has some special properties. This establishes a relation 
between the MINFAS problem and the MINREC problem on Eulerian digraphs, that we explore in the next 
section. 

For two subsets A and B of V, we denote by cutc(A, B) the set {(«, v) 6 E : u e A and v 6 B}. We write 
cutc(A) for cutc(A, V\A), and cut c l (A) for cutc(V\A,A). The following appears stronger than the property 
Vv e V, deg^(v) = deg^iv), but are actually equivalent 

Lemma 1. For every A QV we have \cutc(A)\ = \cut c l (A)\. 

Proof. Let X = {(u,v) e E : v e A},Y = {(u,v) e E : u e A),Z = {(u,v) e E : u e A and v 6 A}. 
We have X = cut^iA) U Z and Y = cut c (A) U Z. Since cut c (A),cut c l (A) and Z are pairwise disjoint, 
\X\ = \cur G \A)\ + \Z\ and \Y\ = \cut c (A)\ + |Z|. Since G is Eulerian, we have = £ (deg c (v) - deg+(v)) = 

veA 

|Z| - \Y\ = \cur G \A)\ - \cut G (A)\. a 

Definition 1. Let A be an acyclic arc set and s a vertex of G. Let rc(A, v) denote the subset of all vertices of 
G that are reachable from s by a path in G[A]. The set A\cut c l {rc(A, s)) U cutc{rc(A, s)) is called cut-stretch 
of A at s. We denote this set by Csc{A, s). 

The idea of cut-stretch is to construct a new acyclic arc set, so that it does not contain less arcs than the old 
one. Moreover, the number of vertices, that are reachable from a fixed vertex, increases after performing 
the cut-stretch. For an intuitive illustration of this definition let us give here an example. Figure la shows 
an Eulerian digraph with an acyclic arc set A shown in Figure lb (plain arcs). If we want to compute the 
cut-stretch of A at V4, we look at all vertices reachable from V4 in G[A]. These vertices are the set rg(A, V4) 
drawn in black on Figure lc. The plain arcs in Figure Id form the set cMf c 1 (rc(A, V4)): arcs of A going from 
the outside (the set {V2, V3, V7}) to rc(A, V4); and the other dotted arcs in this figure form the set cutc{rc(A, V4)) 
: arcs of G going from rc(A, V4) to the outside. Remove the plain arcs in A from A and add the dotted arcs 
of Figure Id, we obtain Csq(A, V4) that is shown in Figure le. 

A simple observation from the above example is that a cut-stretch is still an acyclic arc set and its number 
of arcs is not less than the number of arcs of the old one. The following shows that this property holds not 
only for this example but also holds for the general case. 
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(a) An Eulerian digraph 




(d) the arcs of G going into R from the outside (e) The cut-stretch Csg(A, V4) 

(plain arcs) and the arcs of G going from R to 
the outside (the dotted arcs) 



Figure 1 : An example of cut-stretch 

Lemma 2. Let A be an acyclic arc set and s a vertex of G. Then Csc(A, s) is also an acyclic arc set of G. 
Moreover \A\ < \Cs c (A, s)\. 

Proof. By the definition of cut-stretch there is no arc in Csc{A, s) from a vertex in V\rc(A, s) to a vertex in 
rc(A, s). It implies that if Csc(A, s) contains a cycle, the vertices in this cycle must be completely contained 
either in rc(A, s) or in V\rc(A, s). In this case the arcs of the cycle are also the arcs of A, therefore the cycle 
is also a cycle of A, a contradiction to the acyclicity of A. 

To prove |A| < \Csq(A, s)\, we observe that A n cutc(rc(A, s)) = (from the maximality of rc(A, s)). 
From Lemma 1, we have \Csq(A, s)\ > |A| + \cutc(rc(A, s))\ - \cut c l (rc(A, s))\ = |A|, which completes the 
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proof. 



□ 



The following is the main result of this subsection. 

Theorem 1. Let N be the maximum number of arcs of an acyclic arc set ofG. For every vertex s of G there 
is an acyclic arc set ofN arcs such that it contains no arc whose head is s. 

Proof. Let X be an acyclic set of G of N arcs. We construct a sequence {A,} i£ n as follows: Aq — X and 
A ; = Csc(Ai-i, s) for every i > 1. Lemma 2 and the maximum of N imply that \A,\ = N for every ; 6 N. If 
rciAt, s) = V for some k, A/, is an acyclic set that has the required property since for any vertex v + s of G 
the existence of a path in A* from s to v implies that (v, s) £ A*. The proof is completed by showing that 
there always exists such a k. 

Since a path from s in G[A,] is also a path from s in G[A,+i], we have rc(A,-, s) c rc(A,+i, s). It suffices 
to show that if rc(A,-, s) c V then rc(A;, s) c rc(A, + i, s). Since rc(A,-, s) c V, there is an arc e = (vi, V2) of 
G such that V] 6 rc(Aj, s) and V2 £ ^(A,, s). Since e € A,-+i, there is a path in A,-+i that is from 5 to V2 going 
through vi. It implies that vo e rc(A/+i, 5), therefore rc(A;, s) c rc(A, + i, 5). □ 

Definition 2. A vertex 5 of G is called sink in an acyclic arc set A of G if s has indegree in G[A]. Let A be 
an acyclic arc set of G such that A has exactly one sink s. A vertex s' of G distinct from 5 is called sinkable 
in A if there is an arc of G whose head is s and whose tail is in rc{A, s'). 

We call such a vertex s' sinkable because the idea is to use the arc from s' to s to construct an acyclic arc set 
where it becomes a sink. The fact that this is done by the cut-stretch at s' is stated in the following lemma. 

Lemma 3. Let A be an acyclic arc set ofG having exactly one sink s. If s' is sinkable in A then Csc(A,s') 
has exactly one sink s' and s is sinkable in Csc(A, s'). Moreover A C Csc(Csc(A, s'), s). 

Proof. Let X denote r<j(A, s') and Y = V\X. Since A has exactly one sink s, for any v e V there is a path in 
A from s to v, therefore from s to s'. The acyclicity of G[A] implies that s &Y. 

Clearly, s' is a sink of Cs c {A, s'). To prove that Cs c (A, s') has a unique sink, it suffices to show that for 
any v e V there is a path in Csc(A, s') from s' to v. It is trivial if v € X. We consider the case v e Y. Let 
(v', s) be an arc of G such that v' e X. Such an arc exists because of the assumption of the lemma. By the 
definition of cut-stretch we have (v\ s) e Csq{A, s'). Let Pi and P2 be paths in A from s' to v' and from 
s to v, respectively. It follows from the definition of cut-stretch that Csc(A, s') contains P\. Since v e Y, 
the path P2 goes through only the vertices in Y. Therefore Csc(A, s') also contains Pi. Hence the path 
P\ U {(v', s)} U P 2 is a path in Cs c (A, s') from s' to v. 

Let P3 be a path in A from s to 5'. The acyclicity of G[A] implies that P3 goes through only the vertices 
in Y U {«'}. Therefore there is an arc (v", s') such that v" e F. Clearly, we have rc(Csc(A, s'), s) = T. By 
the definition of the sinkability we have s is sinkable in Cs c {A, s'). 

It remains to show that A c Csg{Csq{A, s'), s). This follows immediately from the fact that 
Cs c (Cs c (A, s'), s) = A U cMf G (T, X). □ 

For each s e V, let Xs denote the number of maximum acyclic arc sets of G with exactly one sink s. 
It is well-known that for an undirected graph G, Tq{\,0) counts the number of acyclic orientations with a 
unique fixed sink, therefore counts Xs, where Tc(x,y) is the Tutte polynomial of G. This implies that if G 
is an undirected graph, Xs is independent of the choice of s. The following is a generalization of this fact to 
Eulerian digraphs 

Proposition 1. For any two vertices s\, s 2 ofG we have x Sl — Xs 2 - 
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Proof. We claim that if (v', v) e E(G) then^,, < Xv- Let 3K\ denote the set of maximum acyclic arc sets of 
G having exactly one sink v, and M2 the set of maximum acyclic arc sets having exactly one sink V . Since 
(v', v) 6 E{G), v' is sinkable in every acyclic arc set in It follows from Theorem 1 and Lemma 3 that the 
the map 6 : 3i\ — > M2, defined by A — > Csc{A, v'), is well-defined. Let A be arbitrary in J{\ . It follows from 
Lemma 3 that A c Csc(Csc(A, v'), v). Since A is maximum, we have A = Csc(Csc(A, v'), v). This implies 
that 9 is injective. Therefore < fflil, equivalently Xv ^Xv- 

The claim implies that for any two vertices v' and v of G such that there is a path in G from v' to v, we 
have ^ v < Xv ■ Since G is strongly connected, there is a path in G from s\ to 52 an d a path in G from 52 to s\ . 
Hence ^ = ^ S2 . □ 

Note that in an undirected graph a maximal acyclic arc set is also a maximum acyclic arc set (and vice 
versa). This fact no longer holds for Eulerian digraphs. The assertion in Proposition 1 is not correct if we 
replace the maximum acyclic arc sets by the maximal acyclic arc sets. 

We recall the definition of the MINFAS problem 

MINFAS Problem 
Input: A digraph G 

Output: Minimum number of arcs of a feedback arc set of G 

When the problem is restricted to Eulerian digraphs, we call it EMINFAS problem for short. Although 
the EMINFAS problem was known to be NP-hard for its multigraph version (Pli 11), it is worth studying 
the computational complexity of the EMINFAS problem since most variants of the MINFAS problem are 
restrictions of the class of digraphs (simple) (see (HMSSY12)). It does not seem that the construction in the 
proof of (Flil 1) is applicable to the case of simple digraphs. By using Theorem 1 and a stronger construction 
we show that the EMINFAS is NP-hard. We work with a general digraph G = (V, E), and construct an 
Eulerian digraph G' so that an optimum value of the EMINFAS problem on G' implies an optimum value 
of the MINFAS problem on G. The graph G' = (V, E') used in the reduction of MINFAS to EMINFAS is 
constructed as follows. 

The basic idea to construct an Eulerian graph G' from G would be to create a new vertex and add arcs 
from this new vertex to any vertex that has more out-degree than in-degree, and arcs from vertices which 
have in-degree greater than out-degree to the new vertex. To avoid multi-graphs, we furthermore add for each 
of those arcs a new vertex in between, which has in-degree and out-degree 1. More precisely, the vertices 
of G are denoted by Vi, V2, ■ • ■ , v„ for some n. If G is already an Eulerian digraph then G' := G. Otherwise 
let G' be a copy of G. We add to G' a new vertex s. For each vertex v, such that deg c {v,) < deg^ivi) we 
add pi new vertices ivy, wy, • ■ ■ , vv, iP . to G', and for each j e [l..pi] we add two arcs (s, vv/j) and (Wjj, v,) 
to G', where pi = deg^{vi) - deg^(vi). For each vertex v, such that deg^{vi) < deg^vi) we add qi new 
vertices w,;i, w,^, ■ • • , w, i9i to G', and for each j e [l..q,] we add two arcs (wy, s) and (v;, Wy) to G', where 
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(a) A digraph G (b) Eulerian digraph G' 




(c) An acyclic arc set of G of maximum cardi- (d) An acyclic arc set of G' 

nality 



Figure 2: Maximum acyclic arc sets 

qi = deg c {vj) - deg c {vi). Formally, the vertex set and the arc set of G' are defined by 
V := {s}UV U y { WiJ :l<j< \deg- G {vd - deg + c ( Vi )\} 

1 <i<n 

E' : = ELI (J {(s, w u ) : 1 < j < deg^yd ~ deg G ( Vi )}U 
|J {(w Uj , vd : 1 < j < deg c { Vi ) - deg G (vd}L> 

deg- G (Vi)<deg+ G (Vi) 

y {(w u , s):l<j< deg~ G (yd - deg G (v t )}U 

deg G (Vj)<deg G (Vj) 

y {(Vi, w u ) : 1 < j < deg c (vi) - deg c (vi)} 

deg G (Vi)<deg G (Vi) 

Figure 2 shows an example of G (Fig. 2a) and the corresponding Eulerian digraph G' (Fig. 2b). Figure 
2c shows an acyclic arc set of G of maximum cardinality. In order to construct an acyclic arc set of G\ we 
add the arcs (s, Wij), (wij, v,) (all the arcs created to offset vertices having out-degree greater than in-degree 
in G) and (v,-, w,j) (half of the arcs created to offset vertices having in-degree greater than out-degree in G) 
to this set, which indeed results in an acyclic arc set of G' of maximum cardinality. The following shows 
that we can always obtain an acyclic arc set of G' of maximum cardinality with this construction. 
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Lemma 4. Let r be the maximum number of arcs of an acyclic arc set ofG, and 

d = 2^ (degc( v i) ~ deg c (\'i)). 

degQ(Vi)<deg+(Vi) 

The maximum number of arcs of an acyclic arc set of G' is 3d + r. 

Proof. The lemma clearly holds if G is an Eulerian digraph, in which case G' = G. We assume otherwise. 
Note that Ad arcs and 2d + 1 vertices are added to G in order to construct G'. Let r' be the maximum number 
of arcs of an acyclic arc set of G' . 

First, we show that 3d+r < r' . Let A be an acyclic arc set of G of r arcs. Let A' = AU{(s, wy) : (s, Wjj) e 
E'} U {(wy, v,) : (Wij, vi) € E'} U {(v/, Wij) : (v;, Wij) 6 £'}. Since A is an acyclic arc set of G and A' contains 
no arc (w;j, s) of A' is an acyclic arc set of G' . The sets {(s, w,j) : (s, w;,/) e £'},{(w !(7 -, v,) : (w;j, v,) e E'} 
and {(v,-, : (v;, wy) e £'} are pairwise-disjoint, and each of them has exactly d arcs, therefore we have 
constructed an acyclic arc set A' of size |A'| = 3d + r. It implies that 3d + r < r' . 

It remains to show that r' < 3d + r. Let B be an acyclic arc set of G' of r' arcs. By Theorem 1 there is an 
acyclic arc set B' of G' of r' arcs such that B' contains no arc (w,j, s) of The set B' must contain all arcs e 
of G' of the form (s, w,- ,•), (w,- ,-, v,) or (v,-, w, ,) since if otherwise, B' U {e} is an acyclic arc set of G' containing 
r' + l arcs. Let A" denote B'\(\(s, w;j) : (s, Wjj) e £'}U{(w,j, v,) : (w, i; ', v,) e E'}U{(v,-, vvy) : (v ; , Wy) e £'}). 
The set A" is an acyclic arc set of G, therefore |A"| < r. It implies r' = = 3c/ + |A"| < 3c/ + r. □ 

A direct consequence of Lemma 4 is a NP-hardness proof for the EMINFAS problem. 

Theorem 2. The EMINFAS problem is NP-hard. 

Proof. Given a general digraph G, the Eulerian digraph G' can be constructed in polynomial time. Let b be 
the minimum number of arcs of a feedback arc set of G', that is, the solution of EMINFAS on G'. Clearly 
\E'\-b is the maximum number of arcs of an acyclic arc set of G'. By Lemma 4 the maximum number of arcs 
of an acyclic arc set of G is \E'\ - b - 3d, where d is defined as in Lemma 4 and is computable in polynomial 
time. Thus the minimum number of arcs of a feedback arc set of G is \E\-(\E'\-b-3d) = b+3d+\E\-\E'\. This 
implies a polynomial-time reduction from the MINFAS problem to the EMINFAS problem. The MINFAS 
problem is NP-hard, so is the EMINFAS problem. □ 



3 NP-hardness of minimum recurrent configuration problem 

3.1 Chip-firing game 

3.1.1 Chip-firing game on digraphs 

Let G = (V, E) be a digraph. A vertex s is called a global sink if deg^(s) = and for any v G V there is a 
path from v to s (possibly a path of length 0). Clearly if G has a global sink then it is unique. 

A configuration c of G is a map from V to N. The value c(v) can be regarded as the number of chips 
stored at v. A vertex v of G is active if c(v) > deg^(v) > 1 . Configuration c is stable if c has no active vertex. 
Firing at v results in the map c' : V — > Z that is defined by 



c'(w) = 



c(w) - deg^(w) ifw = v 

c(w) +1 if v ^ w and (v, w) e E 

c(w) otherwise 



This firing is often denoted by c — » c' . Clearly if v is active then c' is also a configuration of G. In this case 
the firing c — > c' is called legal. If d is obtained from c by a sequence of legal firings (possibly a sequence 
of length 0), we write c — > d. 

A game begining with initial configuration cq and playing with legal firings is called a Chip-firing game. 
Note that at each step of firing there are possibly more than one active vertex, therefore there are possibly 
more than one choice of legal firing. As a consequence, it may be a complicated problem if one wants to 
know the termination of the game. Hopefully, it is not the case for the Chip-firing model since the termination 
has a good characterization. 

Lemma 5. (BL92) Let G be a digraph and c an initial configuration. Then the game either plays forever 
or arrives at a unique stable configuration. Moreover if G has a global sink, the game arrives at a stable 
configuration. We denote by c° this stable configuration. 

3.1.2 Recurrent configuration 

Let G = (V, E) be a digraph with global sink s. Since s is always not active no matter how many chips it has, 
it makes sense to define a configuration on G to be a map from V\\s} to N. In a firing when a chip goes into 
s, it vanishes. Therefore the total number of chips is no longer an invariant under firings. A configuration c 
is accessible if for any configuration d there is a configuration d' such that (d + d') — > c, where d + d' is the 
configuration given by (d + d')(v) = d(v) + d'(v) for any v e V\{s}. Configuration c is recurrent if it is both 
stable and accessible. We denote by REC(G) the set of all recurrent configurations of G. 

Fix a linear order vi < vi < ■ ■ ■ < v„ on V, where n = \V\. The Laplacian matrix A of G with respect to 
the order is given by 

^(d G (vi,Vj) if/*;' 
' J \-deg^(vi) ifi = j 

With the order a configuration can be represented by a vector of Z" -1 , therefore can be regarded as an element 
of the group (Z" _1 , +). Let A\ s denote the matrix A in which the row and the column corresponding to s have 
been deleted. We define an equivalence relation ~ on the set of all configurations of G by c\ ~ ci iff there is 
a row vector z e Z"" 1 such that c\ — c% = z • A\ s . The following shows a relation between the set of recurrent 
configurations and the equivalence classes. 

Lemma 6. (HLMPPW08) The set of all recurrent configurations REC(G) is an Abelian group with the 
addition defined by cffic' := (c + c')°. Moreover, each equivalence class according to ~ contains exactly one 
recurrent configuration, and |REC(G)| is equal to the number of the equivalence classes. 

Naturally, one asks if it is possible to verify efficiently whether a given configuration is recurrent? The 
definition of recurrent configuration does not imply an efficient algorithm for this computational problem. 
Nevertheless, the following implies a polynomial-time algorithm for this problem. 

Lemma 7. (HLMPPW08) Let 6 be the configuration defined by 6(v) = 2deg^(v) for every v e V\{s), and 
e be the configuration given by e(v) — 6(v) — 5°(y) for every v e V\{s}. The configuration e belongs to the 
equivalence class of the identity element, and a configuration c is recurrent iff c = (c + e)°. 

Note that the assertion of Lemma 7 still holds if we replace the definition of 5 in the lemma by 5{v) = deg^.(v) 
for every v € V\{s}. The following is a generalization of Lemma 7, where denotes the zero-configuration, 
i.e. 0(v) = for every v € V\\s} (0 is in the equivalence class of the identity, but is not a recurrent 
configuration). 
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(a) A digraph with global sink s 



(b) A configuration c 




(e)(c + e)° (f)(c+/3)° 
Figure 3: Verifying a recurrent configuration 



Lemma 8. Let Abe a subset of V\[s] satisfying that for every v e V there is a path in G from a vertex in A 
to v. Let jibe a configuration such that ji is in the same equivalence class as and j3(v) > Ofor every v e A. 
Then a configuration c is recurrent iff c — (c + f5)°. 

Proof. =>: Let c — (c+ f3)°. The proof is completed by showing that c is recurrent. Configuration c is stable, 
therefore it remains to prove that c is accessible. Let d be a configuration. Since c is recurrent, there is a 
configuration d" such that c - (d + d")°, therefore c = (c + fj)° - (d + d" + f3)° . Let d' = d" + fJ. We have 
c = (d + d')°. 

<=: For k e N let k/3 be the configuration defined by (k/3)(v) = k-/3(v) for every v € V\{s }. Since for every 
v € V there is a path from a vertex in A to v, with k large enough and by an appropriate sequence of legal 
firings the configuration kfi arrives at a configuration c' that satisfies c'(v) > deg^(v) for every v e 
We have c = (c + k/3)° — (c + c')° . Since (c + c'){v) > deg^(v) for every v e V\{s}, (c + c')° is accessible, 
so is c. □ 

Lemma 7 is a special case of Lemma 8 with A = V\{s} and P - e. 

Figure 3a shows a digraph with global sink s and a configuration c on the right (Fig. 3b). If we want 
to decide whether this configuration is recurrent, we construct the configuration e (Figure 3c), and compute 
the stable configuration (c + e)° (Figure 3e). Lemma 7 states that c is recurrent if and only if c — (c + e)°. 
However e has a large number of chips, and the computation of (c + e)° may be long. It may be more time 
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(a) An Eulerian digraph 



(b) A recurrent configuration 




(c)c+jS 



(d) Firing graph corresponding to the sequence 

(V 5 ,V1,V2,V4,V3) 



Figure 4: An example of firing graph 



efficient to us another configuration with fewer chips, for example the one given on Figure 3d, /3, which 
allows to decide in a similar way if c is recurrent (Lemma 8 applies since the digraph has a global sink s). 
Consider performing the stabilization by hand: one would clearly prefer using f3 to e. 

3.1.3 Chip-firing game on Eulerian digraphs with sink and firing graph 

Let G = (V, E) be an Eulerian digraph (connected) and a distinguished vertex s of G that is called sink. Let 
G\ s be the graph G in which the out-going arcs of s have been deleted. Clearly G\ s has a global sink s. The 
Chip-firing game on G with sink s is the ordinary Chip-firing game that is defined on the graph G\ s . 

Let /? be the configuration defined by for every v e V\M, /3(v) — 1 if (s, v) £ E and f3(v) = otherwise. 
Since G is Eulerian,/? ~ (after firing -1 time every vertex, except the sink). Lemma 8 implies the burning 
algorithm. 

Lemma 9. (Dha90) Configuration c is recurrent if and only if c — (c + /3)°. Moreover if c is recurrent 
then each vertex of G except for the sink fires exactly once during any sequence of legal firings to reach the 
stabilization of(c + f3). 

Note that the configuration c + (3 can be regarded as the configuration resulting from firing the sink in the 
configuration c. Lemma 9 allows to define the notion of firing graph that is originally from (SchlO). 

W\ H>2 W'3 Wk 

Definition 3. Let c be a recurrent configuration and c + J3 — do — > di -4 d% — > d^ ■ ■ ■ — » dk a legal firing 
sequence of c such that di = c. This sequence of legal firings can be presented by {w\,wi, • • ■ , wu) since dj 
is completely defined by w\,Wz, • ■ • , w; for i > 1. Lemma 9 implies that k = \V\ - 1 and {w\,W2, • • ■ , Wk) = 
V\{s). The graph T = («V, S) with <V = V and 8 = {(s, w,) : (s, w,) e E) U {(w,-, wf) : i < j and (w f , w,) e £} 
is called a firing graph of c. 

Figure 4a presents an Eulerian digraph with the sink s in black. Figure 4b presents a recurrent configu- 
ration. The configuration c + /? is presented in Figure 4c. Starting with the configuration c + j3 we can fire 
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Figure 5: A maximal acyclic arc set 



consecutively the vertices V5, vi, V2, V4, V3 of V in this order to reach again c. With the legal firing sequence 
(V5, V\,V2, V4, V3) we have the firing graph that is presented by the undotted arcs in Figure 4d. Note that legal 
firing sequences of c + {! are possibly not unique, so are firing graphs of c. In the next part we are going to 
study a kind of recurrent configurations that always have a unique firing graph. 

3.2 Minimal recurrent configurations and maximal acyclic arc sets 

In this subsection we work with the Chip-firing game on an Eulerian digraph G = (V,E) with sink s. For 
two configurations c' and c we write c' < c if c'(v) < c(v) for every v 6 A recurrent configuration 

c is minimal if whenever c' + c and c' < c, c' is not recurrent. When c has the minimum total number of 
chips over all recurrent configurations, we say that c is minimum. Let M be the set of all minimal recurrent 
configurations of the game. 

Let J{ be the set of all maximal acyclic arc sets A of G such that s is a unique sink of A. Note that 
maximal acyclic arc set can be considered as a generalization of acyclic orientation on undirected graphs. 
Figure 5 shows such a maximal acyclic arc set of the Eulerian digraph shown in Figure 4a. This subsection 
is devoted to showing that if a recurrent configuration c is minimal, c has a unique firing graph and the set 
of arcs of this firing graph is a maximal acyclic arc set. This gives a map from M to M. Moreover we show 
that this map is a one-to-one correspondence between M and £R. The correspondence can be generalized 
easily to the case when G has multi-arcs. 

When G is an undirected graph, the correspondence is exactly the one that was given in (SchlO). The 
correspondence in (SchlO) deals with the case when G has many sinks. However, the many-sink case is 
not harder than the single-sink case since we can contract many sinks to a single sink, and consider the 
contracted graph. This subsection mainly focuses on showing a relation between M and J[, and not all 
results presented here are needed for the proof of the NP-hardness exposed in the next subsection. The 
following shows a basic relation between acyclic arc sets and recurrent configurations. 

Lemma 10. Let A be an acyclic arc set such that s is a unique vertex of indegree in G[A] and A contains 
all vertices of G. Then the configuration c defined by c(v) — deg^ix') — deg c ^(v) for every v 6 V\{s} is 
recurrent. 

Proof. Since G[A] is acyclic, there is a linear order vq < Vj < V2, • • • < v\v\-\ on V such that if (v,-, vj) e A 
then ; < j. Clearly vo = s. The proof is completed by showing that (vi, V2, ••• , V|v(G)|-i) is a legal firing 
sequence of c + /?. Since Vj is an out-neighbor of s in G[A], we have c(vi) = degg{v\) - 1, therefore it is 
active in c + f3. Now by induction, suppose that (vi, V2, • • • , V;) is a legal firing sequence of c + f3, where 
j < |V(G)| - 1. By firing consecutively the vertices Vj, V2, ■ • ■ , Vj ; in this order we arrive at the configuration 
c' . It suffices to show that V/+i is active in c' . It is clear that vy+i receives 2 ^g(v;, v,+i) chips from its in- 

0<i<j 
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neighbors after all vertices Vj , V2, • • ■ , vj have been fired. Since 2 <fc(vi> - deg c ^(vj + \), the number 

0</<; 

of chips stored at vy+j in c' is not less than deg^ivj+x), therefore Vj+i is active in c' . The claim follows. □ 

From the definition of firing graph, a recurrent configuration may have many firing graphs. However, 
the following implies that the numbers of arcs of those firing graphs have a lower bound that depends on the 
recurrent configuration. 

Lemma 11. If c is a recurrent configuration ofG\ s then for every firing graph T — CV, £) of c, s is a unique 
vertex of in-degree and £ is an acyclic arc set of G. Moreover, f is connected and for each v e V\{s} we 
have c(v) > deg^(v) - deg r (v). 

Proof. It follows immediately from the definition of firing graph that s is a vertex of in-degree in T and 
£ is an acyclic arc set. We show that there is no other vertex of in-degree in T ■ Let (vj , V2, ■ • ■ , V|y|-i) be 
a legal firing sequence of c + /3 that is used to construct T ■ By convention vo = s. For each 1 < i < \V\ — I 
let c' denote the configuration obtained from c +/? by firing consecutively the vertices vi, V2, ••• , v,_i. Since 
v, is not active in c but active in c', v,- must receive some chips during this firing process. This implies that 
there is j < i such that v,) e E. It follows from the definition of firing graph that (vj, v,-) € f, therefore 
degyr(\'i) > 1 . Since T is acyclic and has exactly one vertex of in-degree 0, T is connected. 

It remains to prove that for every v e V\{s} we have c(v) > deg^(v) - deg^-(v). For every 1 < i < \V\ — 1 
vertex v, receives deg r (\'i) chips from its in-neighbors when all vertices vi, Vz, ■ ■ ■ , v,_i have been fired. At 
this point v,- is active, therefore c(v,) > deg^(\>j) - deg r (vj). a 

The notion of firing graph gives a map from M to J{ that is shown in the following. 

Lemma 12. Let c e M and f — ("V, S) a firing graph of c. Then c(v) = deg^(v) — deg^iy) for every 
v e V\{s] and & e Moreover, the configuration c contains \E\ — deg^,(s) — \S\ chips. 

Proof. Let c' be the configuration defined by c'(v) = deg^{v) - deg^{v) for every V\{s}. By Lemma 10 c' 
is a recurrent configuration. It follows from Lemma 1 1 that c' < c. Since c is minimal, we have c' = c, 
therefore c(v) = deg^{v) - deg^(v) for every v e V\{.s}. 

To prove 8 e M, we assume otherwise that there is A e J\ such that & c A (from Lemma 1 1 we 
know that & is an acyclic arc set, hence it is not maximal). Let c" be the configuration defined by c"(v) = 
deg]j(v) - deg cw (v) for every v € V\{s}. Let («,«') 6 A\&. Clearly deg c ^(u') > degj-(u'), therefore 
c"(u') < c(u'). It implies that c" + c and c" < c, a contradiction to the fact that c e At. 

The number of chips c contains is 2 c ( y ) = 2 (deg^{v) - deg^(v)) = 2 deg^iv) - deg^(s) - |£| = 

v&s vts veV 

\E\ - degUs) - |£|. The second statement follows. □ 

For two non-repeated sequences f = (vi, Vz, • ■ • , V|v|-i) an d 9 - (wi,m>2, ■ ■ ■ , W|v|-i) of the vertices in 
y\{s), pre£(f, g) denotes the maximum integer k such that for every i satisfying 1 < i < £, we have = Wk- 
Note that if Vi + w\ then pref(f, g) = 0. The following shows that there is a well-defined and injective map 
from M to J{. 

Lemma 13. For every c e M, c has exactly one firing graph. 

Proof. Let fi = (vi, V2, ■ ■ • , V|y|_i) and \z = (w\, Wz, ■■■ , W|v|-i) be two different legal firing sequences of 
c + p. Let j denote pref(fi, f 2 ) and f = (vi, Vz, • ■ ■ , Vj, Wj + \,v j+ \, v' j+v v' j+4 , ■■■ , the sequence of 

vertices of G, where (V +3 , v^ +4 , ■ ■ • , v jV|_i) is the sequence (v^+2, ■ • ■ , V|v|-i) with Wj+\ deleted. Clearly, f is 
also a legal firing sequence of c + /3. Let Ti = (*Vi , £1) and T' - CV' , £') denote the firing graphs of c with 
respect to fi and f, respectively. 
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We claim that T\ — T' ■ Lemma 12 implies that \&\ \ — |£'| = £ deg + c (v)- £ c(v). Hence it suffices 

veV\(s) reV\{i) 

to prove that £i\£' = 0. We assume otherwise that £i\£' + 0. Let k denote the integer such that Wj+\ = \'k- 
Note that k > j+ 1. It follows from the definition of firing graph that £i\£' = {(v,, v*) e E : j+1 < i < k— 1}. 
Let X = {(v,-, V/t) : (v,-, v^) e T'} and Y = {(v,-, v&) : (v,-, v^) e T\\. Since f can be viewed as fi in which \\ 
has been moved backward, we have X Q Y. It follows from £i\£' + that X c Y, therefore deg r ,( v k) < 
deg r (vi), a contradiction to the assertion of Lemma 12. 

Let Ti denote the firing graph of c constructed by f 2 . The proof is completed by showing that T\ -Ti- 
Let 8 = (81,82, ■ ■ ■ , <5|v|-i) be a legal firing sequence of c+fj such that the firing graph constructed by 8 is the 
same as T\ and pref(<5, f2) is maximum. We are going to show that 8 = f2- Let p denote pref((5, f2). If 8 + \% 
then p < \V\ — 1. Let 6' denote the sequence (Si, 8%, • ■ ■ ,8 P , w p+ \,8p + \,u p+ ^, u p+ \, ■ ■ ■ , u\v\-i) of vertices of 
G, where {u p+ 3, u p+ 4, ■ ■ ■ , M|v|-i) is the sequence (8 p +2, 8 p+ 3, ■ ■ ■ , <5|v|— l) with the vertex w p+ \ deleted. The 
above claim implies that the firing graph of c constructed by 8' is the same as the one constructed by 8. It is 
clear that pre£(5', f2) > pref(5, f2), a contradiction to the maximum of pref((5, f2) □ 

For two non-repeated sequences f = (vi, \>2, • ■ • , v\v\-i), 9 = (w\,W2, ■•• , W|v|-i) of vertices in V\{s] we 
denote by inter(f, g) the sequence (vi, V2, • • ■ , Vj, Wk+\, v^+i, v' k+v v' k+4 , ■ ■ ■ , v' m i ), where k = pref(f, g) and 
( v it+3' v 'k+v ' ' ' ' v [v|-i) ^ s tne sec l uence (Vjfc+2) V£+3, ■ • ■ , V|v|-i) with the vertex w^+i deleted. It is easy to see 
that pref(f, g) < pre£(inter(f, g), g). Note that if f and g are two legal firing sequences of a configuration 
c, inter(f, g) is also a legal firing sequence of c. The following result is the converse of Lemma 12. 

Lemma 14. Let A e J{ and f denote G\A\. Then the configuration c defined by c(v) — deg^(v) — deg^iy) 
for every v 6 V\{s} is a minimal recurrent configuration. 

Proof. For a contradiction we assume otherwise that c is not minimal. There is c' e M such that c' + c and 
c' < c. Let T' be the firing graph of c' . By Lemma 12 we have E(f') e J?l and T' + T ^ 

Since A is acyclic, there is a non-repeated sequence fi = (vi, v%, ■ ■ ■ , \'|v|-i) of vertices in V\{s) such that 
if (vi, vj) € A then i < j. Clearly, f i is a legal firing sequence of c + fi. Similarly, there is a non-repeated 
sequence f2 = (w\, W2, • ■ • , W|v|-i) of vertices V\{s) such that if (w,-, wj) e E{T') then ; < j. Clearly, f2 is a 
legal firing sequence of c' + j3. We define the sequence {g,};<=N as follows 

9o = fi 

= inter(g ; , f 2 ), i > 

Let p be the minimum integer such that g p - f 2 . Note that for every / > p, g,- = f 2 . Since T + T' , there is 
a minimum integer q < p such that the firing graph constructed by g^ = (81,82,- ■■ , 8\v\-\) is distinct from 
the firing graph constructed by Q q +i. Let k = pre£(g 9 , f 2 ) and / be the integer such that 81 = w^+i- The firing 
graphs constructed by g„ and g ?+1 are denoted by Q\ and Q%, respectively. 

We claim that for every k + 1 < ; < / - 1 we have (8j, Si) t E. For a contradiction we assume otherwise. 
By a similar argument as in the proof of Lemma 13, the set of arcs of Q2 whose head 8/ is a subset of 
the set of arcs of Q\ whose head 8\. The assumption implies that there is an arc e e E such that e e Q\ 
and e $. Q2, therefore deggJiSj) < degg^Sfi. Since pref(g,, f2) < pre£(g, + i, f2) for every < i < p - 1, 
degg {Si) is equal to the in-degree of 81 in the firing graph constructed by g ; , = f2, namely T' ■ It follows that 
deg r {8i) = deggi^Si) > degg(8{) = deg' r ,(8i), therefore c{8[) < c'{8i), a contradiction to the fact that c' < c. 

Since E(gi)\E{g 2 ) = {(8],8i) e E : k+l < i < l-l], it follows from the above claim that E(gi)\E(g 2 ) = 
0, therefore E{@\) c E(@2)- The choice of q implies that E{@\) — A, a contradiction to the fact that A is a 
maximal acyclic arc set. □ 

The following is the main result of this subsection. 
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Theorem 3. Let f c denote the firing graph of c, the map from Mto J[ defined by c i— > T c is bijective. 



Proof. Lemma 12 and Lemma 13 imply that the map is well-defined and injective. Lemma 14 implies the 
surjectivity. □ 

We end this subsection with an interesting property of the Chip-firing game on Eulerian digraphs 

Proposition 2. The number of minimum recurrent configurations is independent of the choice of sink. 

Proof. Theorem 3 and Lemma 12 imply that the map c i-» T c induces a map from the minimum recurrent 
configurations to the maximum acyclic arc sets of G in JK. Therefore the number of minimum recurrent 
configurations is equal to the number of maximum acyclic arc sets of G in M. It follows from Proposition 1 
that the number of maximum acyclic arc sets of G in 3i is independent of the choice of sink, so is the number 
of minimum recurrent configurations. □ 

Proposition 2 states that the number of minimum recurrent configurations is characteristic of the digraph 
itself. 



3.3 NP-hardness of minimum recurrent configuration problem 

In this subsection we study the computational complexity of the following problem 



MINREC problem 

Input: A graph G with a global sink. 

Output: Minimum total number of chips of a recurrent configuration of G. 



If the input graphs are restricted to undirected graphs G with a sink s, the problem can be solved in poly- 
nomial time since all minimal recurrent configurations have the same total number of chips, namely 
Nevertheless, the problem is NP-hard for general digraphs. In particular, we show that the problem is NP- 
hard when the input graphs are restricted to Eulerian digraphs. 



EMINREC problem 

Input: An Eulerian digraph G with a sink s. 

Output: Minimum total number of chips of a recurrent configuration of G. 



Theorem 4. The EMINREC problem is NP-hard, so is the MINREC problem. 

Proof. Let G be an Eulerian digraph with sink s. Let k be the maximum number of arcs of a feedback arc 
set of G and k' be the minimum number of chips of a recurrent configuration of G. Since the EMINFAS 

problem is NP-hard, the proof is completed by showing that k + k' — £ degt,{v). 

veV\{s) 

By Theorem 1 there is an acyclic arc set A of G such that \A\ = k and s is a unique vertex of indegree in 
G[A]. Lemma 10 implies that the configuration c defined by c(v) = deg^(v) - deg c ^(v) for every v € V\{s} 
is recurrent. Clearly k + £ c ( v ) - 2 degtiy) and k + k' < £ degt(v) since G is Eulerian. 

veV\{s) veV\M veV\{s) 

It remains to prove that k+ k' > £ deg^,{v). Let c be a recurrent configuration such that 2 c{v) - k' . 

veV\\s) veV\{s) 

Let T be a firing graph of c. Lemma 1 1 implies that c(v) > deg + c {v) - deg r (v) for every v e V\{s}, therefore 
k + k'> 2 c{v) + \E{T)\> 2 deg + c (v). □ 

veV\{s] veV\{s) 
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4 Conclusion and perspectives 



In this paper we pointed out a close relation between the MINFAS problem and the MINREC problem. The 
important consequence of this relation is the NP-hardness of the MINREC problem. It would be interesting 
to investigate classes of graphs that are situated strictly between the class of undirected graphs and the class 
of Eulerian digraphs, for which the MINFAS and MINREC problems are solvable in polynomial time. We 
discuss here about such a class. 

It follows from Theorem 1 that to compute the maximum number of arcs of an acyclic arc set of an 
Eulerian digraph, we can restrict to the acyclic arc sets that satisfy the condition in Theorem 1 . With different 
choices of s we have different sets of maximal acyclic arc sets. One would prefer to choose a vertex s 
such that all maximal acyclic arc set have the same number of arcs since a maximal acyclic arc set can 
be computed quickly, therefore a maximum acyclic arc set. Figure 6a shows an Eulerian digraph. If vi is 
chosen, we have exactly one maximal acyclic arc set that is shown in Figure 6b. If V2 is chosen, we have 
exactly two maximal acyclic arc sets with different sizes. Thus one computes easily a maximum acyclic arc 
set if vi is chosen. 




(a) An Eulerian digraph (b) A maximal acyclic arc set with respect to 

vi 




(c) Maximal acyclic arc sets with respect to i>2 



Figure 6: Maximal acyclic arc sets with different choices of s 

Note that there are many Eulerian digraphs in each of which there is no vertex s that satisfies this good 
property. By an experimental observation we see that the class of Eulerian digraphs, for which at least one 
vertex s has the property, is rather large. However, a characterization for this class of graphs, on which 
the MINFAS problem is polynomial, is unknown and remains to be done. In addition, the observation also 
provides a heuristic algorithm for the EMINFAS problem. It is interesting to investigate the properties of 
this algorithm. 

We also presented in this paper a number of interesting properties of feedback arc sets and recurrent 
configurations of the Chip-firing game on Eulerian digraphs. One of the most interesting properties is the one 
in Proposition 2. We propose here an open question that is currently in our interests for further investigations: 
Is there any stronger result for Proposition 2 on Eulerian digraphs, and on digraphs? We believe that the 
results we presented in this paper can be generalized to general digraphs. 
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